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Abstract. 

XMM-Newton background maps for the 3 EPIC instruments in their different instrument /mode/filter combina- 
tions and in several energy bands have been constructed using a superposition of 72 pointed observations. Event 
datasets, with point sources excised, for the different instrument /mode/filter combinations have also been assem- 
bled, with longer exposure times than previously available files. The construction of the background maps and 
event files, together with their properties and usage are described here. Also given are statistical properties of the 
photon and particle components of the XMM-Newton EPIC background, based on the analysis of the 72 datasets. 
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1. Introduction 

The XMM-Newton observatory (Jansen et al. 2001) pro- 
vides unrivalled capabilities for detecting low surface 
brightness emission features from extended and diffuse 
galactic and extragalactic sources, by virtue of the large 
field of view of the X-ray telescopes with the EPIC MOS 
(Turner et al. 2001) and pn (Striider et al. 2001) cameras 
at the foci, and the high throughput yielded by the heav- 
ily nested telescope mirrors. The satellite has the largest 
collecting area of any imaging X-ray telescope. 

In order to exploit the excellent EPIC data from ex- 
tended objects, the EPIC background, known now to be 
higher than estimated pre-launch, needs to be understood 
thoroughly. With a good model of the particle and pho- 
ton background, one can correctly background-subtract 
images and spectra extracted over different energy bands 
and from different areas of the detectors. 

Here we provide details of a project to use a large 
number of XMM-Newton pointed observations to help de- 
fine the EPIC background and to produce background 
maps for each of the three EPIC instruments (pn, MOSl 
& M0S2) in several different energy bands. Also, sig- 
nificantly improved background event files, useful for 
spectral analysis, with longer exposures than previously 
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produced files, and specific to several particular instru- 
ment/mode/filter combinations, have been created as part 
of the analysis. The co-addition of many fields allows the 
minimization of any 'cosmic variance', resulting from vari- 
ations in the local diffuse X-ray emission or contamination 
from pathologically bright sources. 

This paper is intended as an aid when studying or 
working with the EPIC background, and in particular, 
when using these files. The current understanding of the 
XMM-Newton background is described briefly in Section 
2. Section 3 describes the methods used in the creation of 
the background products. Information on how to use the 
background images for XMM EPIC background analysis 
is given in Section 4, along with caveats as to their use. 
The relevant information for the event files can be found 
in Section 5, followed by some concluding remarks. 

All the background product files (maps, event files, re- 
lated software etc.), together with other scripts and pro- 
cedures for XMM-Newton Background Analysis are avail- 
able from http : //www . sr . bham . ac . uk/xmmS/ 

2. The XMM-Newton EPIC X-ray Background 

The EPIC background has been shown (via the work of 
Lumb et al. (2002) and many others; see Appendix) to 
comprise solar soft protons, cosmic rays, electronic noise 
and cosmic X-ray photon background. We now briefly dis- 
cuss the properties of these contributions: 

— Solar soft protons (see Sect. 3. 2) perhaps accelerated by 
'magnetospheric reconnection' events and trapped by 
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the Earth's magnetosphere, then gathered by XMM- 
Newton's grazing mirrors. These dominate the times 
of high background. 

During quiescent periods (i.e. with no significant soft 
proton contamination), the remaining components are; 

— High energy, non- vignetted cosmic ray induced events, 
unrejected by the the on-board electronics. Also, asso- 
ciated instrumental fluorescence, due to the interaction 
of high-energy particles with the detector. 

— Non-vignetted electronic noise: bad (bright) pixels, 
and dark current, though the latter is thought negligi- 
ble. In actuality, most of the bright pixels are rejected 
on board the satellite, and the vast majority of the re- 
maining bad pixels are removed by analysis software 
on the ground. The amount of 'flickering' activity of 
these remaining pixels is thought to be related to the 
radiation level of the observation. 

— Low to medium energy, vignetted X-ray photons from 
the sky. These can be divided into the local (predom- 
inately soft) X-ray background, the cosmic (harder) 
X-ray background, and single reflections entering the 
telecope from bright sources outside of the nominal 
fleld of view (FOV). Lumb et al. (2002) estimate that 
the contribution of difi'use flux gathered from out-of- 
fleld angles of 0.4—1.4 degrees is of order 7% of the 
true, focussed in-fleld signal, and the associated sys- 
tematic error (due mainly to the energy dependence) 
is ±2%. 

Table 1 gives a brief summary of the temporal, spatial 
and spectral properties of these, the major components 
contributing to the XMM-Newton EPIC background. 

The present analysis is related in several respects to the 
work of Lumb et al. (2002), and the reader is encouraged 
to consult this work. Additional notes regarding other re- 
lated work on the EPIC background can be found in the 
Appendix. 

3. Data Analysis 

Source-subtracted, high-background-scrcencd and 
exposure- and source-corrected images (maps) of the 
particle and photon components of the EPIC background 
have been created separately for each EPIC instrument, 
and in several different energy bands. This has been per- 
formed separately per observation, over a large number 
of individual observations. 

The individual background maps for a particular in- 
strument and in a particular energy band have then been 
collected together (for the same instrument mode/filter 
combination) over the whole set of observations. Via vari- 
ous cleaning, flltering and 'cr-clipping' techniques, a 'mean' 
background map is created (for each particular back- 
ground component /instrument/energy band/mode/filter 
combination) . 

Procedures have been written to perform the differ- 
ent aspects of the analysis, making extensive use of the 



Table 2. The coordinate systems of the 1' background 
detector maps and the 4" detector maps used in the anal- 
ysis. 

Instr. Coord. Range in No. of pixels 

(DETX/Y) Detector coords. in maps 







min. 


max. 


(4") 


(1') 


PN 


DETX 


-19199 


14400 


420 


28 




DETY 


-16799 


15600 


405 


27 


Ml 


DETX 


-20399 


20400 


510 


34 




DETY 


-20399 


20400 


510 


34 


M2 


DETX 


-20399 


20400 


510 


34 




DETY 


-20399 


20400 


510 


34 



XMM-SAS tasks, Chandra CIAO tools and HEASOFT's 
FTOOLS utilities. 

Before discussing the analysis in depth, we describe 
some of the files and terms used, and the structure of 
some of the final products. 

3.1. Data Analysis: Files & Definitions 

We first define a number of terms used below: 

— Background maps. Detector maps of all the back- 
ground components combined i.e. soft and hard 
X-ray photons (vignetted), internal particles (non- 
vignetted), some small component of soft protons 
(scattered/funnelled) and single-reflections from out- 
of-FOV sources. 

— Detector maps. Images in detector coordinates (DETX 
DETY) (as opposed to sky coordinates), defined, in the 
present analysis, such that the individual pixels are 1' 
by 1', DETX/DETY of (0,0) lies at a pixel intersection, 
and the full area of the CCDs is covered. A scale of 1 
arcminute has been chosen so as to give us decent count 
statistics in each pixel over the full energy band. More 
detailed information as to the coordinate system used 
(and the 4" detector maps used within the analysis) 
is given in Table 2. Software has been written (see 
Sect. 4.2) to rebin the detector maps to any spatial 
scale, and to reproject the maps to any sky coordinates 
(specified via a user-input sky image). 

— Instruments. The present analysis has been performed 
for each of the EPIC instruments; pn, MOSl & M0S2 
(hereafter PN, Ml, M2). 

— Instrument mode & filter. Several different data acqui- 
sition modes exist for each of the EPIC instruments. 
Also, observations have been taken using different fil- 
ters for each of the EPIC instruments. The present 
analysis has been performed for the most common in- 
strument mode/filter combinations, and those most 
useful to the analysis of diffuse X-ray emission from 
extended objects. These are full-frame mode with thin 
filter (hereafter ft) and full-frame mode with medium 
filter (fm), for each of PN, Ml & M2, and also (for PN) 
full- frame-extended mode with thin filter (et) and full- 
frame-extended mode with medium filter (em). 
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Table 1. Summary of the components within the XMM-Newton EPIC Background; temporal, spatial and spectral 
properties 





PARTICLES 




PHOTONS 




SOFT PROTONS 


INTERNAL 
(Cosmic-ray induced) 


ELECTRONIC 
NOISE 


HARD X-RAYS 


SOFT X-RAYS 


Source 


Few 100 keV 
solar protons 


Interaction of High 

Energy particles 
with detector 


1) Bright pixels 

2) Elcc. overshoot 
near pn readout 


X-ray 

background 
(AGN etc) 


Local Bubble 

Galactic Disk 
Galactic Halo 


Variable? 
(per Obs) 
(Obs to 
Obs) 


Flares (>1000%) 
Unpredictable. 
More far from 

apogee. 

Low-E flares turn 
on before high-E 


±10% 
±10% 

No increase after 
solar flares 


±10% 

1) >1000% (pbcels 
come and go, also 
meteor damage) 


Constant 
Constant 


Constant 
Variation with 
RA/Dec (±35%) 


Spatial 

Vignetted? 

Structure? 


Yes (scattered) 

Perhaps, 

unpredictable 


No 

Yes. Detector + 

construction 

MOS: outer CCDs more 

Al, CCD edges more Si 

PN: Central hole in 

high-E lines (~8keV) 


No 
Yes 

1) Individual pixels 
& columns 

2) Near pn 
readout (CAMEX) 


Yes 
No 


Yes 

No, apart from 
real astronomical 
objects 


Spectral 


Variable 
Unpredictable 
No correlation 
between intensity 
+ shape 

Low-E flares turn 
on before high-E 


Flat ± fluorescence + 
detector noise 
MOS: 1.5 keV Al-K 
1.7 keV Si-K 
det.noise<0.5 keV. 
High-E — low-intensity 
Unes (Cr, Mn, Fe-K, Au) 
PN: 1.5 keV Al-K 
no Si (self-absorbed) 
Cu-Ni-Zn-K (~8keV) 
det. noise < 0.3 keV 


1) low-E (<300eV), 
tail may reach 
higher-E 

2) low-E (<300eV) 


~1.4 power law. 
Below 5keV, 
dominates over 
internal 
component 


Thermal with 
<,lkeV emission 
lines 



— Energy bands. The analysis has been performed 
in the following XMM-XID & PPS standard en- 
ergy bands; energy bands 1—5, defined as follows; 
1: 200-500eV 2: 500-2000eV 3: 2000-4500eV 4: 
4500-7500 eV 5: 7500- 12000 eV, and a full energy 
band 0: 200-12000 eV. 

— Closed datasets. A few XMM observations have been 
obtained with the filter wheel in the 'closed' po- 
sition. No photons reach the CCDs, so the event 
files contain only the instrumental and particle com- 
ponents of the background. Such datasets have 
been collected together and processed, and exist 
on ftp : //www-station. ias .u-psud. fr/pub/epic/ 
Closed (Marty et al. 2002). They are dependent 
on the instrumental mode, and have been processed 
for Ml full-frame (exposure 110655s), M2 full-frame 
(101336 s), PN full-frame (42079 s) and PN full-frame- 
extended mode (26731s). 



3.2. Data Analysis: Preparation and Reduction 

The initial analysis is essentially the same for each obser- 
vation, and is detailed in the following steps: 

— For each observation, the relevant Pipeline Processing 
System (PPS) products (event lists, source lists, atti- 
tude files, housekeeping files etc.), from the standard 
analysis carried out at the SSC, are collected together. 

— For each instrument, region files are created from the 
PPS source lists. These are then used to remove all 
the source events from each of the relevant event files. 
A conservative extraction radius of 36" is used to re- 
move the sources (for comparison, Lumb et al. (2002) 
used 25"). These regions are also removed from previ- 
ously created mask files (these are required to calculate 
losses in area due to source removal). 

— A visual inspection is made of the data to make sure 
that there are no strange features in the field, and 
to ascertain whether there are any very bright point 
sources or large diffuse sources which could contam- 
inate the background, even after source subtraction. 
Datasets which fail this inspection are rejected from 
any subsequent analysis. 
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Fig. 1. The effects of high background flaring: The usable 
fraction of data available after flare-screening is shown for 
all the datasets analysed. 

— Each of the event files are then filtered for periods of 
high background (solar proton flares). Lightcurves are 
created using the whole detector, for single events in 
the energy range 10— 15keV (where no source counts 
are expected, due to the very low effective area of 
the telecope) and with FLAG values as defined by 
#XMMEA_EM or #XMMEA_EP. E.g. for PN, using the fol- 
lowing expression for the XMM-SAS task evselect; 

PI in [10000:15000] &&#XMMEA_EP &&PATTERM==0 
Good Time Interval (GTI) files are created from these 
lightcurves by applying an upper count rate threshold 
of 100 (PN) or 35 (M1/M2) ct/lOOs. The event files are 
then filtered, keeping the low count rate time periods. 
The amount of time lost due to periods of high back- 
ground flaring can be seen in flg. 1, where the distri- 
bution of the fraction of data left after flare-screening 
is shown for all the datasets analysed. 

— The event files are then filtered further. Events with 
energies below 150 eV are discarded. For PN, only sin- 
gles and doubles are retained, for M1/M2, singles, dou- 
bles, triples and quadruples are retained. Finally, the 
event lists are filtered using the #XMMEA_EM/P FLAG 
expressions, excepting that events from outside the 
field of view (out-of-FOV) are kept. 

— In the case of PN, it was necessary to further filter 
the files to eliminate a small number of persistent bad 
(bright) pixels/columns, occurring in many (though 
not all) of the observations. The same pixels were also 
found to contaminate the PN closed datasets. These 
pixels were removed from all the PN datasets. Given 
the large-scale (!') of the final BG maps, this had a 
very small effect, as the loss in area is very small. 



Nevertheless, this difference was corrected for in the 
final maps. The pixels removed were as follows: CCDl 
col.13 & (56,75), CCD2 (46-47,69-72), CCD5 col.ll, 
CCD7 col.34, CCDIO col.61, CCDll (47-48,153-156) 
(50,161—164). No bad pixels were removed from any 
of the MOS datasets (pointed observations or closed 
datasets). The event files are now filtered and have 
had all sources removed. 

— For each of the three instruments, a small-scale (4") 
non-vignetted exposure map (with dimensions given 
in Table 2) is created. From the source-removed mask 
file, an area map (4") is created, containing zero values 
at the positions where sources have been removed, and 
unity values elsewhere. These two maps are multiplied 
and rebinned to form a large-scale (1'; see Table 2) 
'area-times-exposure' map (see Fig. 2). 

— For each instrument (of 3) and each energy range (of 
6), (1') detector maps arc created. Because X-ray pho- 
tons cannot reach the out-of-FOV areas of the detec- 
tor, the events in these areas of the detector are solely 
due to the instrumental and particle components of the 
background. By making use of the Closed datasets, and 
comparing the number of counts in the out-of-FOV re- 
gions of the Closed datasets and the particular obser- 
vation dataset in question, one can separate (for each 
mode-dependent instrument and energy range) the 1' 
detector map into a particle (pa) map, containing the 
particle and instrumental background components and 
a photon (ph) map, containing the photon component 
of the background. This is done, producing 2 maps for 
each instrument and each energy range. 

— The particle images are then exposure-corrected via a 
direct division by the LIVETIME (corrected for peri- 
ods of high background and deadtime, dcadtime being 
times when detector areas are affected by cosmic rays, 
and therefore not available to detect X-rays) , and a di- 
vision by a 1' detector mask of the sensitive area of the 
CCDs (to correct for chip gaps). The photon images 
are exposure-corrected via a division by the appropri- 
ate 'area-times-exposure' map. 

3.3. Data Analysis: Observation, Mode and Filter 
Selection 

The above preparatory analysis has been performed for 
116 XMM-Newton observations, spanning a range in in- 
strumental modes, filters, exposure times, and degree and 
duration of high-background flaring. 

In order to produce the final background maps, the 
observations have been grouped together in terms of in- 
strument/mode/filter etc. We then discarded datasets suf- 
fering from any of the following deficiencies: 

— All observations containing a significantly bright 
source whose wings could still contaminate the back- 
ground after source subtraction. 

— All observations containing a large diffuse source, 
which could contribute to the estimated background. 
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Fig. 2. A non-vignetted (4") exposure map is multiplied by a source-removed (4") mask map and rebinned to create 
a large-scale (!') 'area-times-exposure' map (example is for PN). 



— All observations where the background flaring was of 
such an extent that, after flare-removal, less than 10% 
of the original exposure remained. 

The remaining 72 'clean' observations used in the pro- 
duction of the final background files are summarized in 
Table 3, which lists the revolution number, the RA and 
Dec, the Galactic hydrogen column, and a code giving the 
instrumental mode and filter for, respectively, Ml, M2 and 
PN [f - full-frame mode, e - full-frame-extended mode , t 
- thin filter, m - medium filter]. Also given are a nominal 
Ml exposure time, and the fraction of the exposure time 
remaining after removal of high-background periods and 
accounting for deadtime effects (here the PN fraction is 
given) . 

It is instructive to study some properties of the 72 
clean observations. Fig. 3 shows the distributions of (from 
left to right) the revolution number of the observation, the 
Galactic hydrogen column density, and the live exposure 
time, after cleaning for times of high background. 

Almost all of the observations were taken in a 100- 
revolution period between the end of June 2001 and the 
middle of January 2002. The Galactic hydrogen column 
density distribution shows that the majority of the obser- 
vations are pointed in directions of low to medium col- 
umn. A very small number of observations are pointed 
in directions of very high column. The post-flare-rcmoval 
exposure-weighted mean value of the hydrogen column 
over all the 72 observations is 6.83x 10^° cm~^. The post- 
flare-removal exposure times show a systematically larger 
exposure time for the MOS instruments than for the pn 
(the M2 distribution is essentially identical to the Ml 
distribution shown). This is partly due to the pn being 
slightly more sensitive to proton flaring, but is mainly due 
to the fact that operational overheads generally lead to 
the pn having shorter exposure times (by w5 ks) than the 
MOS. Mean values of the live post-flare-removal exposure 
times are: Ml: 22.0ks (standard deviation 13.8ks), M2: 
22.2ks (s.d. 13.8ks), PN: 14.8ks (s.d. 11.9ks). 



Table 4 summarizes the final cleaned observation in- 
formation in terms of the different combinations of in- 
strument, instrument mode and filter used. The nominal 
exposure time is the sum of the individual LIVETIMEs 
i.e. corrected for periods of high-background and dead- 
time (see however the discussion of the exposure maps in 
Section 5.1). 

3.4. Data Analysis: Cubes and Clipping 

Here we describe the production of a single 'averaged' 
background map from many background maps over dif- 
ferent observations. At any given position in DETX/Y 
coordinates, a subset of the datasets used may be affected 
by contamination or diffuse sources. Such effects can be 
removed by clipping outliers, before taking the mean of 
the remaining 'good' images at this position. 

In order to produce a particular set of background 
maps (in the different energy bands) for a particular in- 
strument/mode/filter combination, the following steps are 
taken. Let us take for example the first entry in table 4, 
i.e. MOSl, full-frame mode, thin filter. 

— For each of the 49 relevant observations, comprising in 
total over 1 Msec of clean, low-background data, the 
0-band (full energy band) 1' photon detector maps of 
the background are stacked together into a 3-D 'image- 
cube' (see Fig. 4) of dimensions DETX, DETY (each 
in the case of MOSl of size 34; see Table 2) and N^j^g 
(of size, in this example, 49). 

— For each DETX/DETY, a statistical analysis of the (in 
this case, 49) values at that particular DETX/DETY 
point is performed. Here, a 'clipcube' is constructed 
(with dimensions identical to the input imagecube), 
containing information (I's and O's) as to which cells 
in the imagecube are within the allowed range and 
which values are not, i.e. which values are 'clipped' (see 
Fig. 5). The allowable range is defined as within some 
number of standard deviations (us) from the mean 
value at that DETX/DETY. For the present analy- 
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Fig. 3. Histogram distributions for the 72 clean observations of (from left to right) the observation revolution number, 
the Galactic hydrogen column, and the live exposure time (after cleaning for high background times) for Ml (thick 
line) and PN (filled region). 



Imagecube 

N=5,6,7... 



N 
i 




DETY 
^ DETX 

Fig. 4. A 3-D 'imagecube' of dimensions DETX, DETY and Ng^g- 



sis, this number of as is set at 1.2. Steps are taken, 
including the prior removal of very strongly aberrant 
pixels, to ensure that the initial mean used to define 
the clip limits is not seriously biased by outliers. The 
mean value at that particular DETX/DETY is then 
calculated from the remaining 'u-clipped' values. 
— Imagecubes are then created, as in the 0-band, from 
the intermediate energy band (bands 1—5) 1' photon 
detector maps of the background. Due to the limited 
statistics in the band 1—5 imagecubes, the broad-band 
clipcube is used to create the cr-clipped band 1—5 im- 
ages. The same procedure is performed for the particle 



images, and for all instrument/mode/filter combina- 
tions. 



4. The Background Maps 

cr-clipped, and averaged, exposure-corrected background 
maps have been created as described above for each in- 
strument/mode/filter combination analysed (of 8) and in 
each energy band (of 6). Photon and particle maps have 
been created separately, as have maps with the two com- 
ponents recombined. This leads to a total of 144 back- 
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Table 3. Summary of the final cleaned and filtered observations used in the production of the EPIC background files. 
Given is the revolution number, the J2000.0 equinox RA and Dec, and the Galactic hydrogen column (in units of 
atoms cm~^) towards these coordinates. Also given is a code giving the instrumental mode and filter for Ml, M2 and 
PN [f - full- frame mode, e - full- frame-extended mode . t - thin filter, m - medium filter], plus the nominal exposure 
time, and the fraction of the pn exposure time remaining after removal of high-background periods. 
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ground maps, and these are shown in Figs 6— 9. The maps 
can be obtained from http : / /www . sr . bham . ac . uk/xmm3/ 
Al_ftOOOO_cphim4Ml .f its is, as an example of the 
file naming, an exposure-corrected photon background 
map. A particle map is named *cpaim*, instead of 
*cphini*, and a particle and photon combined map *cim*. 
The example above is a MOSl map (Ml instead of M2 or 
PN), and is in energy band 4 (of 0—6). The mode is full- 
frame (given by the f ), as opposed to (for PN) extended 
full-frame mode (e), and the filter is thin (t), as opposed 
to medium (m). The six character mode-|- filter string is 
as in Table 3, hence the corresponding PN file is named 



Al_0000ft_cphim4PN.f its. The Al indicates a first re- 
lease of these background maps. 

4.1. The Background Maps: Usage 

When might an observer with their own EPIC data want 
to make use of these background maps? In order to anal- 
yse extended objects, and make background-subtracted 
images, estimate low-surface brightness fiux levels, cre- 
ate radial profiles or perform 1- or 2-D surface brightness 
fitting, one needs a map of the appropriate background. 
Sometimes however, the diffuse, extended nature of the 
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Table 4. Summary of the cleaned and filtered observations used in the production of the EPIC background files, 
separated into the different combinations of instrument, instrument mode and filter used. The exposure time is the 
sum of the individual LIVETIMEs i.e. corrected for periods of high-background and deadtime. The post-flare-removal 
exposure-weighted mean value of the hydrogen column is also tabulated. 
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'Clipcube' 

N=5,6,7... 



N 

A 




DETY 
^ DETX 

Fig. 5. A 3-D 'clipcube' of dimensions DETX, DETY and Ng^g. In this example, for DETX/DETY=(1,1), the 
value in image 1 is outside of the allowable range, and is to be clipped from the calculation of the mean value 
for DETX/DETY=(1,1). 



user's target source is such that the determination of the 
background from their own dataset is difficult or impossi- 
ble. Also, using a background from a significantly spatially 
removed section of the same data leads to problems with 
vignetting and detector variations. Hence the need for the 
independently-produced background maps created here. 

What follows is a recipe for how an observer can make 
use of the background maps in conjunction with their own 
source data. Caveats, and problems that can occur when 
the recipe is not or can not be followed are given there- 
after. 



— The source data should be flare-rejected to a similar 
level to that applied in the creation of the background 
maps. ^ 



— Images should be created from the flare-rejected source 
datasets in the same XID/PPS energy bands (bands 
0—5) as used here. ^ 

— As a particular source image will be created (most 
commonly) in sky coordinates, and will have a reso- 
lution finer than 1', the routine BGrebinimage2SKY 
should be used to rebin and sky-project the equivalent 
background map(s) to the resolution and sky position 
of the source image (see Sect. 4.2). 

— Depending on whether the source maps are raw count 
maps or exposure-corrected flux maps, the equiva- 
lent rebinned background map(s) (which are exposure- 
corrected) may be scaled appropriately, either by an 
exposure time, or using an appropriate exposure map. 

— A comparison of the counts or flux in the out-of- 
FOV areas of the source image and background map 



A.M. Read and T.J. Ponman: The XMM-Newton EPIC Background 



9 






Fig. 6. The final cr-clipped, averaged, and exposure-corrected background maps, for MOSI, full-frame mode. Energy 
band (full-band), then bands 1—5 run from left to right. Rows 1—3 are for medium filter and show (rowl) the 
particle and photon combined (cim) maps, (row 2) the particle (cpaim) maps, and (row 3) the photon (cphim) maps. 
Rows 4—6 show the equivalent maps for thin filter. All 0— band maps are to the same level of scaling. All other maps 
are to a scaling 1/3 that of the 0— band maps. Note that at the very edges of the photon (and combined) map FOVs, 
there are sometimes a few pixels with unusually large or small values. This is due to extremely small exposure values 
in the original observation maps amplifying the noise. 



yields the scaling by which the particle component (the 
'cpaim' map) of the background needs be scaled. 
If source-free regions within the FOV {^xist within the 
source datasct, then count /flux comparisons in these 
regions yield the scaling for the photon component (the 
'cphim' map). The scaled photon component can then 
be add to the scaled particle component to give the 
final background map. ^ 



^ If such a flare-rejection method as used here leads 
to zero or very few Good Times, then the user's data is 
heavily flare-contaminated, and the background maps are 
therefore not suitable for source data extracted from the 
whole dataset. Small discrepancies due to slightly differ- 
ent flare-rejection methods may be accounted for by the 
subsequent scaling. 

^ Though the user should work using the same energy 
bands (bands 0—5) as used in the present analysis, note 
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Fig. 7. The final cr-clipped, averaged, and exposure-corrected background maps, as for Fig. 6, but for M0S2, full-frame 
mode. 



that the standard XID/PPS bands have been used, and 
that extension to larger energy bands can easily be per- 
formed by summing individual band images. For example, 
band 1 is heavily contaminated by detector noise in the 
particle background, so the user may prefer to work in 
band (2-|-3-f 4-1-5). 

^ If source-free regions within the FOV do exist within 
the source dataset, then the user is advised to work as 
above with the particle and photon ('cpaim' & 'cphim') 
maps separately. The particle and photon scaling factors 
may be different (hence the fact that separate particle 
and photon background maps have been made available) . 
If no source- free regions exist, the user may be forced to 
assume that the scalings are the same. Here, just the 'cim' 
images need be used. A script for comparing out-of-FOV 



counts (from event files), compareoutofFOV, is available on 
http://www.sr.bham.ac.uk/xmm3/ Note especially that 
at the very edges of the map FOVs, there are a few pixels 
with unusually large and small values. This is due to ex- 
tremely small exposure values in the original observation 
maps amplifying the noise. These areas should be avoided 
when calculating scaling values. 

4.2. The Background maps: Statistics 

For each of the instrument /mode/filter combinations, and 
for each of the energy bands 0—5, mean count rates (in 
units of ct ks^^ arcmin^^) and standard deviations about 
these mean values have been calculated using all of the 
contributing observations listed in Table 3. This has been 
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Fig. 8. The final cr-clipped, averaged, and exposure-corrected background maps, as for Fig. 6, but for pn, full-frame 
mode. 



done separately for the photons and for particles. For the 
photons, the mean values were calculated within the cen- 
tral 16'xl6', to avoid the small number of pixels discussed 
above with unusually large and small values. Table 5 lists 
these values for the photons, table 6 lists them for the 
particles. 

Several features are evident from the tables. For a par- 
ticular instrument, the photon and particle background 
count rates are of the same order, the values for the pn 
being a factor of a few greater than for the MOS. 

Whereas, for the photons, the decline in intensity (in 
ct ks~"^ arcmin"^ keV~^ with increasing energy is quite 
steep, for the particles, it is far flatter, emphasizing the 
points made in Table 1, i.e. that the particle component 
of the background appears essentially flat, and dominates 



above 5keV over the photon (~ 1.4 power law + soft ther- 
mal lines) component. This can be seen in Figs. 6—9 also; 
at low energies, the photon maps are brighter, whilst at 
higher energies, the particle maps dominate. In studying 
the actual values, one can see that, on average, for the 
MOS, the photons dominate the background over most of 
the spectrum, the particles only attaining a comparable 
count rate in the highest energy band. For the pn how- 
ever, the particles attain a comparable count rate to the 
photons at a lower energy, and are then far more domi- 
nant at the highest energies. Again, this can be seen in 
Figs. 6-9. 

The standard deviation in the mean values represents 
the scatter in the mean background count rates over the 
different observations. Note though there is likely some 
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Fig. 9. The final cr-clipped, averaged, and exposure-corrected background maps, as for Fig. 6, but for pn, extended 
full-frame mode. 



Table 5. Mean count rates for the photon background maps (over the central 16'xl6'), for the different instrument, 
mode and filter combinations, and for each of the six (five plus total) standard energy bands. 



Instr. 


Mode/ 


Nobs 




Mean Count Rate (ct ks ^ 


arcmin ^) (+ standard deviation) 






filter 




Band 0: 


Band 1: 


Band 2: 


Band 3: 


Band 4: 


Band 5: 








200-12000 eV 


200-500 eV 


500-2000 eV 


2000-4500 eV 


4500-7500 eV 


7500-12000 eV 


MOSl 


ft 


49 


2.10 (1.36) 


0.37 (0.18) 


0.78 (0.46) 


0.49 (0.37) 


0.35 (0.22) 


0.31 (0.18) 


MOSl 


fm 


21 


2.00 (1.06) 


0.32 (0.12) 


0.80 (0.38) 


0.48 (0.27) 


0.33 (0.16) 


0.29 (0.14) 


MOS2 


ft 


46 


2.23 (1.40) 


0.39 (0.19) 


0.84 (0.50) 


0.53 (0.39) 


0.37 (0.22) 


0.32 (0.18) 


MOS2 


fm 


26 


1.81 (1.07) 


0.31 (0.12) 


0.76 (0.39) 


0.42 (0.27) 


0.29 (0.16) 


0.27 (0.13) 


PN 


ft 


18 


6.50 (3.89) 


1.90 (1.02) 


2.46 (1.28) 


0.93 (0.71) 


0.69 (0.43) 


0.61 (0.24) 


PN 


fm 


12 


4.31 (2.24) 


1.13 (0.50) 


2.04 (0.94) 


0.72 (0.33) 


0.64 (0.36) 


0.68 (0.48) 


PN 


et 


32 


5.19 (4.38) 


1.87 (2.06) 


1.71 (1.01) 


0.86 (0.68) 


0.72 (0.52) 


0.79 (0.55) 


PN 


em 


8 


5.41 (2.66) 


1.94 (0.77) 


1.90 (0.47) 


0.90 (0.60) 


0.76 (0.50) 


0.74 (0.40) 
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Table 6. Mean count rates for the particle background maps for the different instrument, mode and filter combinations, 
and for each of the six (five plus total) standard energy bands. 



Instr. 


Mode/ 


Nobs 




Mean Count Rate (ct ks ^ 


arcmin ^) (+ standard deviation) 






filter 




Band 0: 


Band 1: 


Band 2: 


Band 3: 


Band 4: 


Band 5: 








200— 12000 eV 


200-500 eV 


500-2000 eV 


2000-4500 eV 


4500-7500 eV 


7500— 12000 eV 


MOSl 


ft 


49 


1.40 (0.11) 


0.12 (0.03) 


0.44 (0.04) 


0.24 (0.02) 


0.26 (0.02) 


0.34 (0.03) 


MOSl 


fm 


21 


1.43 (0.11) 


0.13 (0.03) 


0.45 (0.04) 


0.24 (0.02) 


0.26 (0.02) 


0.34 (0.02) 


M0S2 


ft 


46 


1.34 (0.09) 


0.14 (0.02) 


0.42 (0.03) 


0.23 (0.02) 


0.24 (0.02) 


0.32 (0.02) 


M0S2 


fm 


26 


1.31 (0.09) 


0.13 (0.02) 


0.42 (0.04) 


0.23 (0.02) 


0.24 (0.02) 


0.31 (0.02) 


PN 


ft 


18 


8.37 (2.29) 


2.13 (0.43) 


1.95 (1.36) 


1.50 (0.86) 


1.13 (0.31) 


2.05 (0.21) 


PN 


fm 


12 


8.16 (1.60) 


2.23 (0.31) 


1.55 (0.77) 


1.32 (0.47) 


1.11 (0.20) 


2.08 (0.25) 


PN 


ot 


32 


7.96 (1.52) 


2.10 (0.29) 


1.61 (0.86) 


1.32 (0.58) 


1.15 (0.32) 


2.01 (0.20) 


PN 


em 


8 


8.22 (2.72) 


2.27 (0.87) 


1.58 (0.81) 


1.31 (0.64) 


1.12 (0.31) 


2.05 (0.26) 



contribution from errors introduced in the separating of 

particles from photons, a large contributor being the pres- 
ence of the extremely small exposure values in the orig- 
inal observation maps amplifying the noise in the out- 
of-FOV regions. This effect is more prevalent in the pn 
datasets (where the out-of-FOV regions are smaller). Also, 
there will be some contribution to the scatter due to the 
Poissonian variation in the out-of-FOV counts used to nor- 
malise the particle contribution. For the lower-exposure 
datasets, this can approach 10% in some energy bands. 

As regards the particles, the scatter seen in the MOS 
particle background rates is small, about 10%, and ap- 
pears roughly constant with energy. The pn experiences 
a larger scatter of about 20%, and of this, the greatest 
scatter (50%) is observed at lower energies. Not all of this 
scatter however, is thought due to variations in the actual 
particle count rates (sec above). The photons show quite 
a range in count rate scatter (from 30—60% and beyond), 
with no particular trend with energy or instrument. 

Fig. 10 shows an example of the distribution of back- 
ground count rates observed within the sample. Shown are 
(thick line) the Ml 0-band full-frame thin filter photon 
count rate distribution [49 entries; Table 4] (the equiva- 
lent M2 distribution appears almost identical), and (filled 
region) the PN 0-band full- frame-extended thin filter pho- 
ton count rate distribution [32 entries]. The pn distribu- 
tion has a mean value (in this example) some « 2.5 times 
that for the MOS, and has quite a larger scatter. 

4.3. The Background maps: Sky coordinates, 
Rebinning and Reprojecting 

Usually, of course, the user will work in sky coordinates 
(RA & Dec), not detector coordinates, and to a much finer 
resolution than 1'. With that in mind, software has been 
developed to rebin and reproject onto the sky any pro- 
vided background map (scaled or otherwise) to the spatial 
scale and sky position of a user-input image. 

BGrebinimage2SKY is a shell script plus FORTRAN 
routine to convert the low-resolution DETX/DETY back- 
ground maps into high-resolution sky (X/Y) images by 
interpolation. The user gives a template image contain- 
ing the attitude information and the pixel numbers and 
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Log (Count rate) (Ct/s/arcmin2) 

Fig. 10. (Thick Une) the Ml 0-band full-frame thin filter 
photon count rate distribution, and (filled region) the PN 
0-band full-frame-extended thin filter photon count rate 
distribution 

sizes, and an associated event file (ideally the one used 
to create the template image file) for general header pur- 
poses. A rebinned background map is produced with the 
same resolution and at the same sky position as the tem- 
plate image. The radius of the circle over which interpola- 
tion is performed can be given, (note that the background 
DETX/DETY maps are of 1' resolution, so an interpola- 
tion radius of this or slightly larger is recommended). An 
example of the task in use can be seen in Fig. 11. 

5. Event files 

As part of the analysis, for each observation, event files 
have been created, in the same manner as Lumb et al. 
(2002), filtered for times of high background and with 
all sources removed, as described in Sections. The rele- 
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Fig. 11. Usage of BGrebimmage2SKY to rebin and sky-project a given low-resolution DETX/Y background map (top 
left, specifically the Ml thin-filter 0-band particle background map from Fig. 6) to the spatial binning and sky position 
of a user-created high-resolution sky image (top right). The resultant high-resolution sky background image is shown 
at the bottom. The dark areas in the background map are due mainly to CCD gaps. 



vant event files have been merged together for each instru- 
ment/mode/filter combination (Tables 3 & 4), and eight 
separate event files have therefore been created (with nom- 
inal exposure times given as in Table 4), each having an 
extension containing a calibrated event list in the same 
format as produced by the XMM-SAS. The event files 
have had sky coordinates assigned to them for a pointing 
on the sky of RA=0, Dec=0, PA=0. These event files (or 
indeed any) can be reprojected onto any point in the sky 
via e.g. skycast (see [http://www.sr.bham.ac.uk/xmm3/J . 



It is believed that these event files offer several im- 
provements over previous versions (e.g. those of Lumb et 
al. (2002), themselves improvements on previously-created 
versions) for several reasons. Important points are as fol- 
lows (note that many of these points apply not only to 
the event files, but also to the background maps discussed 
above); 

— The event files have been created separately for each 
combination of instrument, instrument mode and fil- 
ter. Hence, event files now exist for medium filter and 
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also for pn fuU-frame-extended mode, neither of which 
had existed previously. 

— As stated in the analysis sections, all data are collected 
from source-subtracted and high-background filtered 
fields with no bright sources or diffuse features, which 
could contaminate the 'background', even after source- 
subtraction. 

— The datasets are longer than previously available event 
files (over a million seconds of clean low-background 
data exists in each of the thin filter MOS datasets, for 
example). This improves signal-to-noise. 

— The secondary extensions of the files have been re- 
moved, and all headers cleaned and corrected. This 
results in far smaller and more manageable files. It is 
possible to treat the event files as normal SAS event 
files; images, spectra and lightcurves can be created 
via evselect, and one can run ftools on the files. 

— The event files combine a large number (see Table 4) 
of observations (rather than a few, long observations) , 
such that the holes left by the source-subtraction are of 
less importance, and are heavily diluted by data from 
other observations. Due partly to the different expo- 
sure times of the individual observations, no single flux 
limit exists above which sources have been removed. 
However, flux histograms of the detected and there- 
after removed sources do show quite a sharp cutoff at 
PS 1 X 10~^^erg cm~^ s~^) (see Fig. 12). 

— The exposure times are believed to be accurate, 
and incorporate the effects of losses due to dead- 
time. For each event file, 4" resolution exposure 
maps (in DETX/DETY coordinates as described 
in Table 2) are also available from the web URL 
http://www.sr.bhain.ac.uk/xmin3/. These are useful 
in evaluating the degree of the source and bad pixel 
removal, and calculating the exact value of the expo- 
sure of the background event files in a particular region 
of the detector. They are named in a similar manner 
to the event files (e.g. expmap_f tOOOO_Ml . f its is a 
MOSl full-frame, thin filter exposure map). The eight 
exposure maps are shown in Fig. 13. 

5.1. Event files: usage and caveats 

The eight event files are available from 
http://www.sr.bhcmi.ac.uk/xmm3/ and are named 
in a similar manner to the background maps. 
El_ftOOOO_Ml.f its is a MOSl full-frame, thin fil- 
ter event file. PN and M2 are for the pn and M0S2 
instruments, an e denotes extended full-frame mode, and 
an m, medium filter. 

Before using the event files for background analysis, 
it is strongly suggested consulting Lumb et al. (2002). 
Because of the above improvements, a number of the 
caveats noted by Lumb et al. do not apply to the present 
files. Many however, are still valid:- 

— Though the extraction of the background can be done 
in detector coordinates, the background event files can 



be skycast onto the sky position of the user's field (see 
pittpT77www.sr.bham.ac.uk/xmm3/ - alternatively use 
the XMMSAS task attcalc). 

— Remnant low-level fiares will remain in the event files. 
It is known, for instance, that at low energies, the pro- 
ton flares turn on more slowly, but earlier than at 
higher energies. The user could apply further, more 
stringent flare-screening to the event files than that 
described in Sect. 3.2, but note that the lowest level 
proton fluxes may be spectrally variable, so that no 
complete subtraction may be possible. 

— Although point sources have been removed, examina- 
tion of images created from the datasets does reveal 
fluctuations in intensity. Over scales of arcminutes ap- 
propriate to extended sources, this is not expected to 
be a significant problem, and indeed representative of 
unresolved background. If a user however, should try 
to extract spectra from regions comparable to the PSF 
scale, care should be taken, and a manual inspection 
may be necessary to guard against a local excess or 
deficit of counts arising from the point source extrac- 
tion procedures. In some of the present cases, it is true 
that, as many of the observations were pointed such 
that the target source was at the centre of the detec- 
tor, a deficit in counts (by as much as 50% or greater) 
is seen at this position. This can be corrected for by 
making use of the exposure maps provided (Fig. 13). 

— Many defects are seen at the lowest energies (below 
0.3 keV), and the calibration of the EPIC response at 
these energies is not so well understood as it is at 
higher energies. Care should be taken when performing 
analysis at ;<0.25keV. 



6. Future Developments 

At the time of writing, there has been only limited expe- 
rience gained in using the maps, the event files and the re- 
lated software. Any feedback on using the datasets would 
be very useful and is most welcome (please email any com- 
ments to £mir30@star . le . ac. uk) . 

It is hoped that further releases of the datasets, created 
using larger numbers of pointed observations, and perhaps 
using further modes and filters, will be made available in 
the future. This will be announced in the usual manner 
and via the URL http://www.sr.bhajn.ac.uk/xmm3/. 
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Fig. 12. Flux histograms of the detected and thereafter removed sources for each of the instrument /mode/filter 
combinations. (Top) Thin filter: MOSl, M0S2 and pn full-frame, pn full-frame extended. (Bottom) the same but for 
medium filter. All plots are to the same scale. 
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Appendix: The XMM-Newton 
Background: Additional notes 



EPIC X-ray 



There have been various studies of the XMM-Newton 
EPIC background, and many of these are described in 
Marty (2002), Lumb (2002) & especially Lumb et al. 
(2002). 



A number of datasets have also been analysed and 
merged together, in a similar manner to that presented 
here. Two sets of data are particularly relevant to the 
present paper: 

(i) Background datasets have been produced (Lumb 
2002, Lumb et al. 2002) for the three EPIC instru- 
ments by source-subtracting and co-adding a few long 
observations. These event files can be obtained from 
http : //xmm . vilspa . esa . es/ccf /epic/#background . 

(ii) Event lists combining several CLOSED ob- 
servations have been created by Phillipe Marty, 
and these have proven very useful in the analy- 
sis and background-subtraction of extended objects 
(see Sect. 3. 2). These data can be obtained from 
ftp : //www-station. ias .u-psud . f r/pub/epic/Closed 

In addition, several novel methods have been used to 
analyse very extended and diffuse X-ray sources, where the 
extraction of the background is difficult. Many of these are 
described in Marty et al. 2002. Of special interest are the 
works of Arnaud et al. (2001, 2002), de Luca & Molendi 
(2002), Ghizzardi et al. (2000), Molendi (2001), and Pratt 
& Arnaud (2002). 
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Fig. 13. 4" resolution DETX/DETY (sec Table 2) exposure maps for the background event files in each of the 
instrument /mode/ filter combinations. (Top) Thin filter: MOSl, M0S2 and pn full-frame, pn full-frame extended. 
(Bottom) the same but for medium filter. 



